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ABSTRACT 

The 'Diffuse Ionized Medium' (DIM) comprises a significant fraction of 
the mass and ionization requirements of the ISM of the Milky Way and is 
now known to be an energetically significant component in most normal 
star-forming galaxies. Observations of the ionized gas in starburst galaxies 
have revealed the presence of gas with striking similarities to the DIM in 
normal galaxies: relatively low-surface-brightness and strong emission from 
low-ionization forbidden lines like [SII]AA6716,6731. In this paper we analyze 
Ha images and long-slit spectra of samples of normal and starburst galaxies 
to better understand the nature of this diffuse, low-surface-brightness gas. We 
find that in both samples there is a strong inverse correlation between the Ha 
surface-brightness {T^Ha) and the [SII]/Ha line ratio at a given location in the 
galaxy. However, the correlation for the starbursts is offset brightward by an 
order-of-magnitude in Ha surface-brightness at a given line ratio. In contrast. 
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we find that all the galaxies (starburst and normal alike) define a universal 
relation between line ratio and the relative Ha surface brightness (Sj^a/Ee, 
where Eg is the mean Ha surface brightness within the galaxy half-light radius). 
We show that such a universal correlation is a natural outcome of a model in 
which the DIM is photoionized gas that has a characteristic thermal pressure 
(P) that is proportional to the mean rate of star-formation per unit area in the 
galaxy {I^sfr)- Good quantitative agreement with the data follows if we require 
the constant of proportionality to be consistent with the values of P and TIsfr 
in the local disk of the Milky Way. Such a scahng between P and T,sfr may 
arise either because feedback from massive stars heats the ISM or because T,sfr 
is determined (or limited) by the mean gas pressure. 

Subject headings: ISM: structure - ISM: kinematics and dynamics - Galaxies: 
ISM - Galaxies: Starburst 



1. INTRODUCTION 

The 'Diffuse Ionized Medium' (DIM) is recognized as an important component of the 
ISM in galaxies. This gas, first discovered as the 'Reynolds Layer' in our Galaxy (see, 
e.g. Reynolds 1990 for a review), seems to be ubiquitous in late-type spiral (e.g. Hoopes 
et al. 1996; Wang, Heckman & Lehnert 1997-hereafter WHL) and irregular galaxies (e.g. 
Martin 1997). Indeed, the universal existence of the DIM has been inferred from integrated 
emission-line ratios of a large sample of normal late-type galaxies (Lehnert & Heckman 
1994). The observed properties of the DIM are characterized by relatively strong low 
ionization forbidden lines compared to normal HII regions, a low surface brightness, a rough 
spatial correlation with HII regions, and a significant contribution ( ~20%-40%) to the 
global Ha luminosity. 

The observations raise many interesting questions about the physical and dynamical 
state of the DIM that are still to be answered. The energy required to power the DIM 
suggests that the gas cither soaks up nearly 100% of the mechanical energy supplied 
by supernovae and stellar winds, or the topology of the interstellar medium must allow 
roughly 1/3 of the ionizing radiation produced by massive stars to escape HII regions and 
propagate into the disk. Although the diffuse nature of the DIM suggests that it maintains 
pressure balance with the rest of the ISM, httle observational evidence has been collected to 
support this idea. Only in the Reynolds layer has the electron density been derived through 
observations of pulsar dispersion measures (e.g. Reynolds 1993) and seems consistent with 
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the typical ISM thermal pressure of ~3000 K cm~^ (Jenkins et al. 1983). Furthermore, it 
is still unknown how or if this pressure is regulated by the hotter coronal-phase gas created 
by supernovae. 

While these questions wait to be answered, a study of the general properties of the 
DIM in galaxies with widely differing star formation rates (SFRs) per unit area or per unit 
volume might provide us with more clues in this endeavor. The simple reason behind this 
is that star formation has a significant impact on the ISM. The energetics and dynamics 
of the DIM must therefore be strongly influenced (or even regulated) by the feedback 
from star formation (e.g. Elmegreen & Parravano 1994). In starbursts, for example, 
the intense star formation can provide feedback to the ISM that might be significantly 
different qualitatively from the case in quiescent galaxies. Dynamically, the collective effect 
of supernovae exploding in a hot rarefied medium created by previous supernovae may 
minimize radiative losses and thus provide the energy to drive a galactic scale outflow in 
starbursts (e.g. Heckman, Armus & Miley 1990). Thus, we would expect a more significant 
kinematic disturbance in the DIM in starbursts than in normal disks and a greater role of 
shock-heating in the energetics throughout the ISM. That is, the hot coronal-phase gas in 
starbursts should be more pervasive than in normal disks (as parameterized by the porosity, 
McKee & Ostriker 1977) and would therefore have the potential to regulate the pressure of 
the ISM over a relatively larger volume. Since the DIM traces the heating and ionization of 
gas occupying a substantial fraction of the volume of the ISM and comprising much of its 
mass, the above issues are closely tied to our understanding of the DIM. 

We note that among the galaxies that have been searched for a DIM, most are normal 
spirals and irregulars. However, observational data suggest that there is a similar DIM 
component in starburst galaxies. For instance, in the nearby starburst NGC 253, Hoopes 
et al. (1996) found a faint Hct-emitting gas surrounding bright HII regions that is similar 
to the DIM in other quiescent spirals. More generally, the spectra of extra-nuclear regions 
of starbursts (e.g. Lehnert & Heckman 1995) show that the relative strengths of the 
low-ionization lines are high — very similar to what has been observed in the DIM in normal 
galaxies. It is therefore worthwhile to systematically study the faint emission-line gas in 
starbursts and compare this gas to the DIM in normal galaxies. 



2. SAMPLE SELECTION AND DATA DESCRIPTION 

In this paper we will explore the emission-line properties of galaxies with a wide range 
of star formation rates (SFRs) using long-slit spectroscopic data. Our sample includes both 
normal, quiescent late-type spirals and IR-selected starbursts. Our data on normal galaxies 
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come from a survey of the DIM in nearby face-on spirals. This subset involves seven galaxies 
selected on the basis of Hubble type (Sb and later), proximity (closer than 10 Mpc), large 
angular size (> 10 arcmin), and relatively face-on orientation (inclination < 65 degrees). 
Part of the data has been presented in Wang et al. (1997), and the rest can be found in 
Wang (1998). The reader can refer to these references for details of the observations and 
data reduction. Complete observational information is summarized in Table 1. 

The starbursts are a sample of infrared-bright (Seo/xm > 5.4 Jy), edge-on (a/b ^ 2), 
infrared-warm (SgOyum/Sioo/um > 0.4) galaxies compiled by Lchncrt (1993, dissertation, 
Johns Hopkins University) and also described by Lchncrt and Heckman (1995). Analysis of 
both samples of data can help us understand the physical state and dynamics of the ionized 
gas over a wide range of star formation rates. 

We note that due to different sample selection criteria, the starbursts are generally much 
more distant and more inclined than our normal spirals. Otherwise, the major difference 
between the starbursts and the normal spirals is just the much larger star- formation rates 
per unit area (Hsfr) in the former (Lehnert & Heckman 1996b). 

In order to compare our normal galaxies and starbursts, we used three quantities that 
have been measured with long-slit spectroscopy at many locations within each galaxy for 
both groups of galaxies. These are the hne ratios represented by [SII]AA6716,6731/Hq;, 
the Ha surface brightnesses, and the Ha or [NII]A6584 linewidths. The [SII] doublet 
was selected to represent low ionization lines rather than [Nil] because Nitrogen has a 
secondary nucleosynthetic origin, while Sulfur has a primary one. That is, the S/0 ratio is 
independent of metallicity, while N/0 oc metallicity (cf. Vila-Costas & Edmunds 1993). As 
shown by WHL this results in systematic differences in the [Nil] /Ha ratio in the DIM in 
galaxies of different metallicity - differences that are not present in the [SII] /Ha ratio. 

For our normal galaxies, we have multiple slit positions that spcctroscopically sample 
representative regions in the disks of these galaxies (WHL; Wang 1998). We find that the 
emission-line properties in these galaxies vary spatially in a smooth way from the centers of 
HII regions out into the surrounding DIM. 

Data for the starbursts are in the form of long-slit spectra centered on the nuclei and 
oriented along both the minor and major axes. The spectra have been extracted using 3 
pixel spacing starting from the nuclei and moving outward along the slit. We can then 
examine the variations in starbursts from bright starburst cores to fainter surrounding 
nebulae, by analogy with our analysis of the normal galaxies. The starbursts are generally 
much more distant than our normal spirals. Therefore we sample the gas in starbursts on 
a much larger scale, as each spatial resolution element would encompass a large number of 
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HII regions. 

The [NII]A6584 linewidth is used to represent kinematics in starbursts, while either Ha 
hnewidth (if EQW(Ha) > 3A) or [Nil] linewidth (if EQW(Ha) < 3A) is used (see WHL) 
for our normal galaxy sample. This difference in the adopted linewidth is unlikely to affect 
our analysis significantly since the Ha and [NH] linewidths correlate well with one another 
for the starburst (Lehnert & Heckman 1995) sample. 

The number of galaxies we can use in the starburst sample is limited by the availability 
of photometric data and other relevant information. We use 32 galaxies from this sample 
that have been observed spectroscopically under nearly photometric conditions and hence 
can provide us with line ratios, Ha surface brightnesses, and linewidths. We have excluded 
from this analysis the Circinus galaxy, whose Ha image in Lehnert & Heckman (1995) is 
dominated by the central AGN and NGC 5253, which is a dwarf galaxy unlike the members 
of our normal galaxy sample. 

In the following analysis wc will attempt to normalize the star-formation rate in a given 
location of a galaxy by dividing the measured Ha surface brightness at that location by 
the average Ha surface brightness in the galaxy. To do so, we define a galaxy effective Ha 
surface brightness Eg as the ratio of half of the total Ha flux to the solid angle subtended 
by the area within the Ha half-light radius(7rrg). This is effectively the average surface 
brightness within the Ha half-light radius. In order to do this scahng we further selected 
19 objects out of the subsample of 32 galaxies for which Lehnert and Heckman (1995) have 
measured the Ha half-light radii and total Ha fluxes. The measured total Ha-|-[NII] flux is 
corrected to the Ha flux based on the measured [Nil] /Ha ratio within 2re. Thus, Eg can be 
estimated for these 19 starburst galaxies and used to scale the observed values of T,Ha as 
measured at different positions in each galaxy. 

All the spectral data in the normal spiral subsample have been properly calibrated for 
absolute surface brightness. In addition, we have measured from the Ha images (WHL 
and Wang et al. 1998 in preparation). Therefore all data from this subsample can be 
utilized in this analysis. 

We have not attempted to correct the observed Ha surface brightnesses for the effects 
of internal reddening. This is likely to be signiflcant. For example, Kennicutt (1983) 
estimates an Ha extinction of about 1 magnitude for typical giant HII regions in normal 
galaxies, and Armus et al (1989) find a typical value of about 2 magnitudes for IR-sclected 
galaxies. Note that while an extinction-correction would increase the absolute values of 
the Ha surface brightnesses (and implied star-formation rates per unit area) it will only 
affect the normalized surface-brightnesses defined above if the extinction varies spatially. 
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We return to these issues later in the paper. 

3. RESULTS 

The relative strengths of the low-ionization emission-lines in starbursts correlate well 
with Ha surface brightness (Figure 1). Data for HII regions and the DIM in our normal 
galaxies are also plotted for comparison. Both sets of data show higher strengths of 
low-ionization lines at lower surface brightness, and they all suggest a strong continuity in 
physical state between the high surface brightness and low surface brightness gas. In fact, 
the correlations can be roughly described as a power- law relation between [SII]/IIa and 
^Ha with similar slopes for both samples. However, there is a noticeable offset between the 
two groups of galaxies: at a given line-ratio, the gas in the starbursts has an average surface 
brightness that is about an order-of-magnitude higher than the gas in the normal galaxies. 
This is not surprising as the starbursts have much higher Ha surface brightnesses in general, 
and correspondingly higher "Esfr- We will further examine this systematic difference in 
surface brightness in the following paragraphs. We reject some DIM data points near the 
nucleus of M 81 due to heating processes other than pure photoionization by young stars in 
that region (cf. Devereux, Jacoby, & CiarduUo 1995, 1996). 

While the hne-ratios correlate strongly with surface-brightness, the starburst sample 
shows no correlation between linewidth and Tina, and only a weak correlation between 
linewidth and [SIIJ/Ha. Linewidths of the gas in starbursts are typically a few hundred 
kilometers per second, about a factor of 10 larger than those in normal galaxies. The data 
for the normal galaxies bear little resemblance to that of the starbursts with respect to these 
kinematic relations. It is likely that kinematics of the emission-line gas are different in the 
two types of galaxies for two major reasons. For the distant highly-inclined starbursts, the 
spectra sample a large number of emission line nebulae across starburst disks even within 
a small aperture, and the lines may be broadened by relative motions among the nebulae. 
Supernova-driven superwinds are also responsible for much of the linewidth broadening 
(Lehnert and Heckman 1996a), especially along the minor axis. Thus we conclude that 
linewidths do not track well with other emission-line parameters for the combined galaxy 
sample. 

In starbursts the continuity between high surface brightness regions (starburst nuclei) 
and low surface brightness (outer) regions, together with the lack of correlation between 
kinematics and emission line intensity ratios, suggests that the gas is mainly photoionized 
instead of mechanically heated. While we have drawn the same conclusion for the gas in 
the normal galaxies (WHL), Figure 1 shows a more general trend that is applicable to 
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the gas in starbursts as well. It is not surprising that most of line emission results from 
photoionization since the ionizing radiation energy from OB stars is about an order of 
magnitude higher than the total kinetic energy released by supernovae and stellar winds 
(cf. Leitherer & Heckman 1995). 

To further explore the possibility of photoionization, we need to understand the 
Ha surface brightness offset between the two groups of galaxies in Figure 1, which is 
presumably physically related to the much higher star formation rates per unit area in 
the IR-sclcctcd starbursts. To test this idea, we adopt the effective surface brightness Sg 
defined in section 2 above to scale the observed values of T^ho and then plot the line-ratio 
[S11]/Hq; versus relative (dimensionless) surface brightness T^Ha/^e (Figure 2). We have 
excluded the starburst data beyond 2re because those data are relatively noisy and are more 
hkely affected by superwind-driven shock-heating processes in the galaxy halos (Lehnert & 
Heckman 1996a). 

Figure 2 shows remarkably universal pattern of line-ratio variation as a function of 
normalized surface brightness. The scaling by Eg has successfully eliminated the systematic 
difference in surface brightness at a given line-ratio between the two groups of galaxies. 
Figure 2 demonstrates the interesting similarities in emission-line properties not only 
between HII regions in normal spirals and in starburst nuclei, but also between the DIM in 
normal galaxies and the relatively faint emission-line gas surrounding starburst nuclei. The 
universal continuity between high surface brightness and low surface brightness gas suggests 
that these emission-line properties all vary with a single parameter and therefore provides 
further support to the idea of photoionization as the dominant mechanism. 

Because of the similarity in emission-line properties, we propose that the fainter 
emission-line gas in starbursts and the DIM in normal galaxies has the same physical 
nature. While it remains to be confirmed obscrvationally that the low surface brightness gas 
in the starbursts is generically diffuse rather than from the sum of many faint, discrete HII 
regions, our conjecture is supported by Hct imaging results of the nearest starburst galaxies 
hke M 82, NGC 253, NGC 5253, and NGC 1569 (cf. Lehnert & Heckman 1995, Marlowe et 
al 1995; Martin 1997). In addition, we (Wang 1998; Wang et al. 1998) have analyzed our 
Ha images of the normal galaxies and find that there is no absolute Ha surface brightness 
limit that cleanly separates the DIM from bright HII regions. Instead, we show that defining 
the DIM in terms of a normalized surface-brightness analogous to that described above 
leads to a natural segregation that is independent of a galaxy's mean surface brightness in 
Ha. This implies that even in a starburst galaxy a DIM component would exist, but this 
diffuse gas will have a characteristically high absolute surface brightness. Since the mean 
TiHa seems to be a reasonable dividing point between the DIM and HII regions, we suggest 
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that for the sample galaxies discussed in this paper, a crude surface brightness limit is the 
T^Ha averaged within the half-light radius {—Tig). Therefore TiHa/^e = 1-0 could be used to 
isolate the 'DIM' gas in Figure 2. 

For simplicity, in the following discussion we will tentatively adopt the same acronym 
DIM for the faint gas in the starbursts. We proceed to use a generic photoionization model 
to explain the correlation in Figure 2 and then address other implications of the data. 

4. DISCUSSION 

4.1. Photoionization of the DIM 

The observed inverse correlation between the relative strength of the low-ionization 
lines and the Ha surface brightness could have a simple physical explanation if the 
emitting gas clouds all have roughly the same density. In this case, the higher the value 
of the local intensity of the ionizing radiation field, the higher the value of the ionization 
parameter U (defined to be the ratio of the densities of ionizing photons and electrons 
within a photoionized gas cloud). Simple ionization equilibrium arguments show that U 
determines the ionization state of the gas, while recombination means that the Ha surface 
brightness will be proportional to the intensity of the ionizing radiation field. Thus, the 
proportionality between U and the Ha surface brightness will naturally produce enhanced 
relative intensities of low-ionization lines in the faint gas (cf. Domgorgen & Mathis 1994). 

The generalized relation between line-ratio and TiHa/^e for the DIM in Figure 2 could 
then be explained i/ there is a direct proportionality between the average thermal pressure 
in the diffuse interstellar medium and the average star-formation rate per unit area in the 
galaxy. That is, the ratio of the Ha surface brightness at a particular location compared to 
the mean value in the galaxy would then be proportional to the local value of the intensity 
of the ionizing radiation field divided by a quantity that is proportional to the thermal 
pressure and hence the density of the photoionized cloud (T ~ 10"^ K). Thus, Tual^e oc U. 

Later in this section, we will briefly discuss some of the physics that might lead to such 
a proportionality. Here we simply remark that there is good empirical evidence that this 
proportionality is roughly obeyed when extreme starbursts like M 82 are compared to the 
disks of normal spirals hke the Milky Way. In the M 82 star burst, the thermal gas pressure 
is P/k ~10^ K cm~^ (Heckman, Armus, & Miley 1990) and the star-formation rate per 
unit area is T,sfr ~ 30 Mq yr"^ kpc~^ for a Salpeter IMF extending from 0.1 to 100 M© 
(cf. Kennicutt 1998). Heckman, Armus, & Miley (1990) and Lehnert & Heckman (1996b) 
show these are typical values for both parameters in extreme starbursts. In comparison, in 
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the local Milky Way disk the thermal gas pressure is P/A; ~ 10^-^ K cm ^ (cf. Jenkins et 
al 1983; Reynolds 1993), and J^sfb. ~ 4 xlO'^ M© yr'^ kpc^^ (McKee & Williams 1997 
adjusted to our adopted IMF). Thus, T,sfr is roughly 7000 times greater in M 82 and the 
pressure is roughly 3000 times greater. This agrees with Lord et al (1996) who estimate 
that both the thermal pressure and FUV intensity in M 82 are three-to-four orders of 
magnitude higher than in the local ISM. 

Let us for the moment then adopt the conjecture that P oc T^sfr, and derive a relation 
between T^Ha, U, and Ue utilizing photoionization models. Suppose the gas in the DIM is 
illuminated by an isotropic ionizing radiation field. Then the one-sided incident ionizing 
flux ^Lyc is related to the observed area-averaged Ha surface brightness of the cloud by 

(Vogcl ct al. 1995) where ena is the ratio of the intrinsic (extinction-corrected) and observed 
Ha surface brightness and f^a is the fraction of recombinations which produce Ha photons 
(=0.46 for T = 10^ K and Case B recombination). The ratio of observed area to total area 
Aproj/Atot is determined by the nebular geometry. It is 1/2 for a slab and 1/4 for a sphere. 
We use 1/3 to represent an average case. According to the definition of U, we can express 
U — A^Lyc/ueC, therefore 

^Ha = 5.9 X 10"^^ (r^^) ^ ^^^^ ^^'^ arcsec-2 (2) 

To compare this to our data, we need to understand how rig in the DIM is related to 
the mean star- formation-rate per unit area, as measured by Eg. Since photoionized gas 
is generically in thermal equlibrium at a temperature of roughly 10^ K (e.g. Osterbrock 
1989), relating to is equivalent to determing the constant of proportionality in the 
relation P cx T^sfr- To do this, we will adopt a purely empirical approach for the moment 
and insist that this constant agree with values for P and T^sfr in the ISM of the Milky 
Way. Later we will explore the possible physical basis of this. 

In the local disk of the Milky Way, T^sfr ~ 4 xlO~^ M0 yr~^ kpc~^ implies an average 
intrinsic Ha surface brightness for the disk of 1.2 x 10^^^ ergs s"^ cm~^ arcsec"^ (where we 
have assumed contnuous star-formation with a Salpeter IMF extending from 0.1 to 100 Mq 
- Leitherer & Heckman 1995). A thermal pressure of P/A; ~ 10^'^ K cm~^ imphes n^, ~ 0.16 
cm~^ in the photoionized gas. Thus, the predicted relation between Eg and Ue based on our 
own Galaxy is: 

Ee = 7.4 x 10"^^ [ ejjl^ ergs s"^ cm-^ arcsec^^ (3) 
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Photoionization models of the DIM (Sokolowski 1993) are able to reproduce the 
observed emission-line ratios provided that the cosmically- abundant, refractory elements 
(i.e. Fe and Si) are largely locked-up in dust grains, as in the case of diffuse clouds in 
our own Galaxy (cf. Savage & Sembach 1996 and references therein). The models also 
better match the data if the radiation field incident on the DIM has been hardened due to 
radiative transfer en route to the DIM (e.g. there is an optical depth of-order unity at the 
Lyman edge between the DIM and the O stars). Adopting these 'depleted and hardened' 
models, we then estimate an empirical relation 



appropriate for the range of the observed DIM lineratios ([SII]/Hq; ^ 0.3 - 1.5). This 
enables us to relate [SIIJ/Ho; approximately to J^Haf^e- The ratio of equations [2] and [3] 
imply that T^Ha/^e = 80 U. 

Thus, using Equation [4] above, we obtain: 



This relation is represented by the solid line in Figure 2. The data agree reasonably well 
with the prediction in the faint gas (e.g. Y^ual^e < !)• The deviation of the data from the 
prediction at higher surface brightnesses {T,Ha/^e > 1) will be briefly discussed in section 
4.3 below. 

As a 'sanity check' we now estimate the values of n^, in normal disks and starbursts 
that are implied by the measured values of Eg based on Eq. [3] . One should keep in mind 
that Eg may need to be corrected for inclination, so the values given here are upper limits, 
especially for the starburst sample where the inclination is high. We infer an average rig of 
0.4 eHa cm~^ from the observed Eg for the normal galaxies. Taking a typical extinction 
of 1 magnitude for Ha (Kennicutt 1983) we obtain rie = 1.0 cm^'^. This is several times 
larger than the value rig ~ 0.16 cm^'^ for the Reynolds layer (Reynolds 1993). There might 
be two major reasons for this discrepancy. Firstly, the Ha extinction in the DIM may be 
less than the typical HII region value of 1 magnitude. Secondly, simple considerations 
of hydrostatic equilibrium (see below) imply that the total ISM pressure (e.g. the sum 
of thermal, turbulent, cosmic ray, and magnetic pressures) decreases with galactocentric 
distance, so that the thermal pressure and hence electron density in the DIM may be 
higher in the inner regions of galaxies. Now, rig for the Reynolds layer has been measured 
in the solar neighborhood (about 8 kpc from the Galactic center), while in our normal 
galaxies is typically 2-4 kpc (see Table 2). The electron density for the DIM in starbursts 



[5//]AA6716,6731 



= 1.1 X 10-^ U 



-0.58 



(4) 




-0.58 



(5) 
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estimated from Ee averages ~ 24 cm~^, after correcting for two magnitudes of extinction in 
Ha (Armus, Heckman, & Miley 1989). This value for rig is then considerably higher than 
that in normal galaxies (as expected). 

To summarize, we have shown that a model in which the DIM in both starburst and 
normal galaxies is photoionized gas whose thermal pressure is proportional to the mean 

rate of star-formation per unit area in the galaxy can quantitatively reproduce the observed 
unified correlation between the ionization state of the DIM ([SII]AA 6716,6731/Ha line 
ratio) and the relative surface-brightness of the DIM shown in Figure 2. To make this test 
we have fixed the constant of proportionality in the relation P oc T^sfb. to its value in the 
local disk of the Milky Way. We now turn to the possible physical basis of this relation. 



4.2. A Supernova-Regulated ISM Pressure 

Suppose we assume that the average thermal gas pressure P within is maintained 
by the energy and mass released by supernovae (and stellar winds) inside r^. We can then 
relate P and therefore in the DIM to the effective surface brightness Sg. 

Chevalier and Clegg (1985) have shown that for the case of spherical symmetry and 
adiabatic conditions 

P = 0.12 ^-2 ^g-) 

Where M and E are the rates at which gas in the ISM is shocked and heated respectively. 
While this is exact for a spherically-symmetric case, its difference from a disk geometry can 
be shown to be negligible (provided that the gas is adiabatic - see below). 

The starburst models of Leitherer & Heckman (1995) predict a simple scaling between 
the rate at which a starburst would return kinetic energy [E) and ionizing photons (Q). 
Case B recombination gives the scaling from Q to the Ha luminosity. Now M is the amount 
of mass per unit time that is heated by supernovae and stellar winds. This will be larger 
than the ejecta directly returned from the massive stars by a 'mass-loading' factor m. For 
a standard Salpeter IMF extending up to 100 M© and a constant rate of star-formation for 
a time longer than 40 Myr, the Leitherer & Heckman models and equation 6 above then 
imply that rig is related to Eg by 

Ee ~ 1.4 X 10"^^ (r^^) arcsec"^ (7) 

where we have assumed P/2nf.k = 10^ K. Eq. [7] agrees with the scaling relation between 
Eg and Ug (eq. [3]) of the local disk, for an appropriate m value (see below). 
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We can compare equation [6] to conditions in the local disk of our own Galaxy. 

Based on the rates at which stellar winds and supernovae inject mass and kinetic energy 
(~8xl0~^ M© yr^^ kpc~^ and 1.2x10^^ ergs s^^ kpc~^ respectively) within 3 kpc from 
the Sun (Abbott 1982; Jura & Kleinmann 1989), the predicted thermal pressure of the 
hot gas is ~2200 m^'^ K cm ^, compared to the representative value of 10^'^ K cm ^ from 
observations of neutral (Jenkins et al. 1983) and ionized (Reynolds 1993) diffuse gas near 
the Galactic midplane. The required amount of mass heated per unit time is about twice 
as much as that injected directly by supernovae and stellar winds (m ~ 2). Only about 1/3 
of this returned mass comes from high- mass stars (M> 5 M©) ,with the bulk coming from 
intermediate-mass AGB stars (Jura & Kleinmann 1989). The situation in starbursts-where 
the mass is returned almost entirely by high-mass stars (Leitherer & Heckman 1995)-is 
therefore somewhat different. Values for m ~ 3 to 10 have been estimated in starburst 
galaxies based on the mass, luminosity, and temperature of the X-ray emitting gas (e.g. 
Suchkov et al 1996; Delia Ceca et al 1997; Wang et al 1997). 

Equation [6] assumes that the thermal gas pressure is determined by the deposition of 
mass and energy by supernovae and stellar winds, that the hot gas that results permeates 
the region of star-formation, and that radiative losses are negligible. These assumptions 
may be valid in starbursts driving superwinds (cf. Heckman, Lehnert, & Armus 1993), 
but probably not in the ISM in normal galaxy disks where the interaction between stellar 
ejecta and the ISM is more complex (cf. Cioffi & ShuU 1991). We therefore consider next 
a different physical interpretation of the relation between gas pressure and star-formation 
intensity. 

4.3. Hydrostatic Equilibrium and Pressure-Regulated Star Formation 

As discussed above, the results in Figure 2 can be understood if the DIM in both 
normal and starburst galaxies is photoionized, has a roughly constant characteristic density 
in each galaxy, and the characteristic thermal pressure in the DIM is proportional to the 
rate of star-formation per unit area in that galaxy. In section 4.2, we considered the 
possibility that the physical coupling was provided by the energy and mass deposited in the 
ISM by supernovae and massive stars. Here, we consider a different interpretation, namely 
that the total pressure in the ISM is specified by hydrostatic equilibrium (e.g. Boulares & 
Cox 1990), and that the star-formation rate per unit area is related to, or perhaps limited 
by, this pressure (cf. Dopita 1985). 

In a case of simple hydrostatic equilibrium, the total (thermal, turbulent, cosmic ray. 
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plus magnetic) mid-plane pressure in a disk galaxy is given by 

Ptot OC T,gas{^* + ^gas) T,gas^tot (8) 

On empirical grounds, it is well-established that the star-formation-rate per unit area 
C^sfr) in disk galaxies scales with both E^a^ and 'Etot- This suggests that there might be a 
simple, direct scaling between T^sfr and Ptot- In fact, Dopita & Ryder (1994) parameterize 
the problem as 

EspR oc E-,Ej;, (9) 

and find empirically that m + n — 2.0 ± 0.5. Kennicutt (1998) finds that T^sfr oc E^;^^, 
while Figure 1 in Dopita & Ryder (1994) implies T,sfr oc E°-^. Except in the most extreme 
starbursts, it is reasonable to take Titot ^ ^gas, so that this last result means roughly that 
^SFR oc E°J. Combining these results suggests that: 

J:SFR oc Ej„tS^? oc Ptoti^gas/^totf-' (10) 

Since there is only a small observed variation in Ylgas/^tot (factors of a few) in the 
disks of late-type galaxies and typical starbursts, this implies that there should be a 
proportionality between Ptot and T^sfr- If we now assume that the thermal component 
of the pressure scales with the total pressure, this is just what we require in order to 
understand Figure 2. 

As we emphasized in section 4.1 above, the rough quantitative agreement between 
photoionization models and the properties of the DIM in the starbursts and normal galaxies 
is independent of the nature of the physical, causal connection between the thermal pressure 
in the DIM and T^sfr- The agreement shown in Figure 2 is based simply on requiring that 
the constant of proportionality between these two quantities is consistent with values in the 
local disk of our Galaxy. 



4.4. The High Surface-Brightness Gas 

While the model of a photoionized, roughly isobaric DIM provides a satisfactory 
quantitative match to the data on the low surface brightness gas (T^Ha/'^e < 1), Figure 2 
shows that there is a systematic offset between the predictions and the data in the high 
surface brightness range (Sj^Q,/Se ~ 1 - 10). This disagreement is in the sense that the 
observed emission-line ratio [SII]/Hq; is too high, and therefore that the actual ionization 
parameter U in the gas must be smaller than predicted. This could be explained if the 
density in the high-surface-brightness gas is higher than estimated in the model. This 
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is entirely plausible, since the high surface-brightness gas (the HIT regions) will hkely be 
significantly ovcr-prcssured with respect to the surrounding diffuse ISM (e.g. the HII regions 
may be self-gravitating or expanding into the lower-pressure DIM). More quantitatively, 
we note that the offset between the model and data for log([SIl]/HQ;) in the bright gas (a 
difference of ~ 0.6 dex on average) would translate into a difference of a factor of ~10 lower 
U and hence higher ng. 

An additional factor is that the Sokolowski models we have utilized for the DIM: f) 
adopted a dust-dcplctcd abundance pattern and 2) assumed that the ionizing radiation field 
had been hardened as it propagated to the DIM. Partial depletion onto grains may occur 
in Hll regions (cf. Garnett et al 1995), but the assumption of spectral hardening is not 
appropriate for the HII regions. Dropping these assumptions would decrease the predicted 
ratio of [SII]/Hq; for a given U, and make the discrepancy worse in Figure 2. This would 
require a decrease in U (and increase in rig) by an additional factor of ~ 3. 

Using the value for rig in the DIM in normal galaxies estimated from Eg above, we 
would then require a density of ~30 cm"^ in the HII regions. This agrees reasonably 
well with the average values of Ue in disk HII regions of ~10-100 cm~^ (e.g. O'Dell and 
Castaheda 1984; Kennicutt, Keel & Blaha 1989) measured with the [SII] and [OH] doublets. 
Our measurements of [SII]A6716/[SII]A6731 for bright HII regions suggest similar values for 

Following the same reasoning, we would estimate that the required density in the 
high-surface brightness gas in the centers of the starbursts must also be about 30 times 
higher than in the low-surface-brightncss gas: ~ 700 cm""'. This is in satisfactory 
agreement with the directly measured central densities of 300 to 1000 cm~^ (e.g. Heckman, 
Armus and Miley 1990; Lehnert and Heckman 1996a). 

5. CONCLUSIONS 

We have compared the emission-line properties of the low surface brightness gas in 
starburst galaxies with the DIM in normal spirals. Both samples show similar attributes of 
enhanced low-ionization forbidden-line strengths (as represented by [Sll]/H«) relative to 
typical Hll region values and a strong inverse correlation between Ha surface-brightness 
and the [S11]/Hq; line ratio (in the form of a smooth transition from high surface brightness 
to low surface brightness regions). The gas kinematics show no strong correlation with 
surface brightness and line ratio in the combined samples. 

Although the Ha surface brightness corresponding to a given [SII] /Ho; line ratio is 
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preferentially about an order-of-magnitude larger in starbursts than in normal galaxies, we 
have demonstrated that this can be understood as a consequence of the proportionately 
higher mean Ha surface brightnesses of the starbursts. That is, we have shown that the 
relative surface brightness at a particular location, defined as the absolute surface brightness 
there (T^Ha) scaled by the mean surface brightness within the Ha half-light radius (Sg) for 
the galaxy as-a- whole, exhibits a remarkably universal correlation with the [Sll]/Ha line 
ratio for normal and starburst galaxies alike. This suggests that the emission-line properties 
of the low surface brightness gas in both groups of galaxies can be unified to a simple 
relation between line ratio and relative surface brightness, and that the variations in line 
ratio and relative surface-brightness are controlled by a single parameter. 

We have constructed a simple photoionization model to explain the correlation between 
T,Ha/^e and hne ratio. We have pointed out that the [SII]/Hq; line ratio has an inverse 
dependence on the ionization parameter U (the local ratio of ionizing photons and electrons 
in the photoionized gas). For simple recombination, TiHa is proportional to the local 
intensity of the ionizing radiation field, //the average thermal pressure in the diffuse ISM 
in a galaxy (P) is proportional to the average rate of star-formation per unit area {Hsfr), 
then since T^sfr can be measured by Sg, it follows that U oc J^Ha/^e- We have argued 
that this result naturally explains the universal dependence of line ratios on relative surface 
brightness. 

Our simple model is able to quantitatively reproduce the data for normal and starburst 
galaxies provided that the constant of proportionality in the relation between P and T,sfr 
is consistent with the observed values for both quantities in local Galactic ISM. Thus, we 
have emphasized that the agreement between our simple photoionization model and the 
data is independent of the detailed physical connection between star formation and ISM 
pressure. 

We have discussed two ways in which P might be physically related to T^sfr- Following 
Chevaher & Clegg (1985) we have first assumed that P is regulated by the feedback of 
mass and energy from supernovae and massive stars. Scaling the amount of mass heated 
per supernova so that the predicted thermal pressure matches the observed pressure in 
the local Milky Way disk, we found that the photoionization model agrees roughly with 
the DIM data in both starburst and normal galaxies. As an alternative, we explored the 
possibility (e.g. Dopita 1985) that T^sfr is determined (or limited) by the total (thermal, 
turbulent, cosmic ray, plus magnetic) pressure Ptot, and that Ptot and P are determined by 
a simple hydrostatic equilibrium condition in galactic disks, i.e. P cx Ptot oc T^tot^gas- We 
used recent empirical results from Kennicutt (1998) and Dopita & Ryder (1994) to argue 
that TiSFR oc Pfqas-i where fgas is the fractional gas mass in the disk. Since igas varies only 
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by small factors, this means that T,sfr does roughly scale with P. 

The simple model can not account for the emission-hne ratios in the high surface- 
brightness gas (the giant HIT regions) unless the densities and thermal pressures there are 
roughly 30 times larger than in the DIM. We argue that this is both reasonable physically 
and in agreement with measurements. 

We conclude that the low surface brightness gas in the starbursts shares a common 
nature with the DIM in the normal galaxies, and propose that the former can be regarded 
as the same gas phase as the latter. Further morphological observations of the low surface 
brightness gas in starbursts can confirm this suggestion. 

We thank S. Baum, D. Calzetti, R. Kennicutt, R. Wyse, C. Norman, C. Martin, and 
A. Ferguson for useful discussions and an anonymous referee for constructive suggestions. 
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Fig. 1. — The [SII]AA6717,6731/Hq; line ratio vs. the Ha surface brightness at the same 
position in the galaxy. The starburst data are represented by small triangles and include 

32 galaxies for which photometric measurements are available (Lehnert & Heckman 1995). 
These data cover the regions from starburst nuclei to their outer disks. The data for normal 
galaxies are represented by different symbols as follows. Stars — M 101, triangles — M 51, 
solid circles— M 81, crosses— NGC 4395, open circles— NGC 2403, diamonds— NGC 6946, 
T — IC 342. These data are from the DIM, bright Hll regions, and intermediate regions 
around isolated HII regions. The DIM points lie to the upper left, the HII regions to the 
lower right, and the intermediate points lie in between. 



Fig. 2. — The [SII]AA6717,6731/Hq; line ratio vs. the relative Ha surface brightness T^Ha/^e 
where Eg is the mean surface brightness within the Ha half-light radius rg. Symbols are the 
same as Figure 1. The data includes only 19 starburst galaxies (in the Lehnert & Heckman 
1995 sample) that have measured values for Eg. Regions beyond 2re are excluded because 
the measurement errors are large and shock-heating by superwinds may be significant in 
these outer regions. This figure demonstrates the universal relation between line-ratio and 
relative surface brightness of the diffuse emission-line gas for galaxies with different SFRs. 
The solid line is the prediction from a simple photoionization model for the DIM assuming 
the thermal pressure in the diffuse gas is related to the star-formation rate per unit area (see 
text for details) . The assumptions in the photoionization models are not appropriate for the 
high-surface-brightness gas — the HII regions with TiHa/^e > 1- 
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Table 1. High- Resolution Spectroscopic Observations^ 



D • • b 

Position 




Oovcrage 


• d 

Grating 


Resolution 


PA 


Blocking Filter 


otandard Star 
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(A) 


\ ) 








M 51: 


















nuc 


6563 


6180-6950 


KPC-24 


0.8-1.0 


60 


GG495 


BD -f262606 


02/17/96 


nuc 


4850 


4450-5250 


KPC-24 


1.2-1.5 


60 


BG39 


Feige 34 


12/04/96 


92" N 40" W 


6500 


6120-6880 


KPC-24 


0.9-1.1 


60 


GC496 


Feige 34 


01/31/97 


92" N 40" W 


4900 


4500-5300 


KPC-24 


1.0-1.2 


60 


4-96 


HZ 44 


01/31/97 


92" N 40" W 
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KPC-24 


1.0-1.2 


60 


4-96 
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02/03/97 
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KPC-24 


0.8-1.0 


60 
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BD -f262606 


02/18/96 


M 81: 
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02/17/96 
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60 
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12/02/96 


316" S 54" E 


6600 


6220-6980 


KPC-24 


0.8-1.0 


60 


GG495 


Feige 34 


12/03/96 


316" S 54" E 


4850 


4450-5250 


KPC-24 


1.2-1.5 


60 


BG39 


Feige 34 


12/04/96 


234" N 183" W 


6500 


6120-6880 


KPC-24 


0.9-1.1 


60 


OG196 


Feige 34 


01/31/97 


231" N 183" W" 


1900 


1500-5300 


Ki'C-2i 


1.0-1.2 


()() 


l-9() 


(;i91I!2B 


02/03/97 


M 101: 


















nuc 


6563 
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KPC-24 


0.8-1.0 
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BD -1-262606 


02/17/96 
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4900 
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KPC-24 
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140 


4-96 


HZ 44 


01/31/97 


286" S 202" W 


6500 


6120-6880 


KPC-24 


0.8-1.0 


140 


GG495 


HZ 44 


02/03/97 


286" S 202" W 


4900 


4500-5300 


KPC-24 


1.0-1.2 


140 


4-96 


G191B2B 


02/03/97 


NGC 2403: 
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nuc 
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KPC-18C 


1.2-1.5 


35 


GG385+CUS04 


Feige 34 


12/02/96 
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4900 


4500-5300 


KPC-24 


1.0-1.2 


35 


4-96 
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02/03/97 
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KPC-24 
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35 


GG495 


Feige 34 


12/03/96 
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KPC-24 
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35 


BG39 


Feige 34 


12/04/96 
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KPC-24 
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35 


4-96 
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02/03/97 
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KPC-24 


0.9-1.1 


35 


GG496 


Feige 34 


01/31/97 


NGC 4395: 
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BD -1-262606 


02/17/96 
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12/02/96 
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60 
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Feige 34 


12/03/96 
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KPC-24 
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60 
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12/04/96 
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KPC-24 
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60 


GG495 


BD -f262606 


02/18/96 


NGC 6946: 
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KPC-24 


0.8-1.0 


140 


GG495 


Feige 34 


12/03/96 


nuc 


4800 


4360-5230 


KPC-18C 


1.2-1.5 


140 


GG385-I-CUS04 


Feige 34 


12/02/96 


35" S 175" W 


6600 


6220-6980 


KPC-24 


0.8-1.0 


140 


GG495 


BD 284211 


12/04/96 


IC 342: 


















nuc 


6600 


6220-6980 


KPC-24 


0.8-1.0 
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GG495 


Feige 34 


12/03/96 


nuc 


4800 
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KPC-24 


0.8-1.0 
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GG495 


Feige 34 
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150 


GG385-I-CUS04 


Feige 34 


12/02/96 



''All observations were made with the KPNO 4m telescope using the RC spectrograph and the T2KB CCD. A slit of 5' 
long and 2" wide is used centered on either the nucleus or the indicated offsets from the nucleus of each object. 

''Centers of the slit. Numbers represents offsets from the galactic nucleus. 

'^Useful spectral range estimated from 1500 pixels of the CCD. 



''All gratings are set in second order. 
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Table 2a. Effective Ha Surface Briglitness 



Object 












(L0) 


(kpc) 


(Mpc) 


(ergs cm~^ arcsec"^) 


M 51 


4.4x10^ 


4.2 


8.4 


2.4x10-16 


M 81 


1.5x10^ 


4.2 


3.6 


6.0x10-1^ 


M 101 


6.3x10^ 


10.2 


7.4 


7.0x10-1^ 


NGC 2403 


1.4x10^ 


2.1 


3.2 


3.2x10-16 


NGC 4395 


1.9x10^ 


1.7 


2.6 


6.0x10-1^ 


NGC 6946 


1.2x10^ 


4.5 


5.9 


5.6x10-16 


IC 342 


2.3x10^ 


2.9 


2.1 


2.2x10-16 



^Integrated Ho; luminosity corrected for foreground Galactic extinction. Data are from 
WHL and Wang et al. 1998 (in preparation). 

''Half-light radius in Ha measured at an inner circle in the images that enclosed one-half 
of the total Ha hght from the galaxy. 

'^Distance adopted from WHL and Wang et al. 1998 (in preparation). 

'^Mean Ha surface brightness within half-light radius Ve derived from Lna (Col. 2) and rg 
(Col. 3). The values have been corrected for foreground Galactic extinction. The values can 
be converted to emission measure based on the relation that 5xl0-i^ erg s"! cm"^ arcsec"^ 
corresponds to 25 pc cm-6^ assuming an electron temperature of 10^ K. 
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Table 2b. Effective Ha Surface Brightness 



Object 


LHa* 










(L0) 


(kpc) 


(Mpc) 


(ergs s~^ cm~^ arcsec"^) 


NGC1134 


6.9x10'^ 


3.7 


46.0 


6.0x10-1'' 


IR03359+1523 


1.6x10* 


1.4 


141.7 


l.OxlO-i"* 


NGC1511 


7.6x10^ 


0.98 


18.0 


9.6x10-1*^ 


NGC1572 


4.0x10'^ 


3.0 


80.5 


5.4x10-16 


NGC1808 


1.2x10'^ 


0.32 


13.6 


1.4x10-1* 


NGC2966 


i.exio'^ 


0.35 


31.7 


1.6x10-" 


IC564 


5.6x10^ 


5.6 


85.6 


2.1x10-1'^ 


NGC3044 


1.4x10'^ 


1.9 


21.7 


4.5xlO-i'5 


NGC3511 


l.lxlO'' 


2.3 


19.4 


2.5x10-1® 


NGC3593 


4.1x10" 


0.49 


12.9 


2.0x10-15 


UGC6436 


7.6x10^ 


2.6 


143.5 


1.4x10-15 


NGC4433 


6.1x10^ 


2.6 


44.8 


1.1x10-15 


NGC4527 


2.9x10^^ 


5.9 


27.7 


1.0x10-"' 


NGC4666 


7.5x10'' 


5.0 


24.8 


3.6x10-"' 


NGC5073 


1.2x10^" 


3.0 


40.9 


1.7xlO-i'3 


NGC5104 


2.5x10''' 


2.6 


79.2 


4.4x10-16 


NGC5775 


3.7x10^^ 


4.0 


25.2 


2.8x10-1® 


ZW049.057 


4.6x10*5 


0.71 


49.7 


1.1x10-15 


IG5179 


1.0x10* 


3.3 


42.3 


1.1x10-15 



^Integrated Ha luminosity corrected for foreground Galactic extinction. Data are from Ha + [NH] luminosities listed by 
Lehnert & Heckman (1995) and have been corrected to pure Ha luminosities based on average [Nil] /Ha ratios within 2 half-light 
radii. 

''Half-light radius in Ha measured at an inner circle in the images that enclosed one-half of the total Ha light from the galaxy 

(Lehnert & Heckman 1995). 

"^Distances adopted from Lehnert & Heckman 1995. 

"^Mean Ha surface brightness within half-light radius re derived from Lhq (Col. 2) and re (Col. 3). The values have been 
corrected for foreground Galactic extinction. The values can be converted to emission measure based on the relation that 
5x10-1^^ erg s-i cm-^ arcsec-^ corresponds to 25 pc cra~^, assuming an electron temperature of 10* K. 
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